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If advanced civilizations appear in the universe with an ability and desire to expand, the entire 
universe can become saturated with life on a short timescale, even if such expanders appear rarely. 

Our presence in an apparently untouched Milky Way thus constrains the appearance rate of galaxy- 
spanning Kardashev type III (K3) civilizations, if it is assumed that some fraction of K3 civilizations 
will continue their expansion at intergalactic distances. We use this constraint to estimate the 
appearance rate of K3 civilizations for 81 cosmological scenarios by specifying the extent to which 
humanity is a statistical outlier. We find that in nearly all plausible scenarios, the distance to the 
nearest visible K3 is cosmological. In searches for K3 galaxies where the observable range is limited, 
we also find that the most likely detections tend to be expanding civilizations who have entered the 
observable range from farther away. An observation of K3 clusters is thus more likely than isolated 
K3 galaxies. 


I. INTRODUCTION 

It is difficult to imagine a scientific discovery that 
would have a more profound impact than direct observa¬ 
tional evidence of advanced civilizations engaged in engi¬ 
neering at the scale of entire galaxies - the so-called Kar¬ 
dashev type III (K3) civilizations [133 ■ Not only would 
such an observation answer age-old questions about intel¬ 
ligent life, but it could also become a source of new infor¬ 
mation about the limits of technology and physics mi- 
Nevertheless, this version of SETI has only recently be¬ 
gun to attract serious attempts at observation mmmm 
125] . with the largest of these searches to date [5] sensi¬ 
tive to technology-induced waste heat from ss 10 5 nearby 
large and dwarf galaxies, and reporting a null result. 

Kardashev’s scale m was advanced in the 1960’s with 
the hope of informing searches of extraterrestrial life, in¬ 
cluding searches for galaxy-spanning civilizations. Now, 
more than 50 years later, there remain very few quan¬ 
titative tools to estimate the number of K3 civilizations 
that could be within range of observation. The problem 
is particularly acute now, with search results beginning 
to be reported - if n nearby galaxies have been examined 
for K3 civilizations with null result, what is the inter¬ 
pretation? Should we have expected to see anything in 
n galaxies to begin with? Based on what assumptions? 
Here we introduce a hypothesis that, if valid, seems to 
heavily constrain the range of possibilities, allowing us to 
make quantitative predictions on the most likely type of 
positive search result: 

Expansion Hypothesis: K3 civilizations 
have, by definition, already exhibited the nec¬ 
essary technology and behavior characteris¬ 
tics required to expand rapidly beyond the 
boundaries of their home galaxy, and are thus 
probable to do so. 
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By “probable,” we mean a probability that is not sup¬ 
pressed by many orders of magnitude. For example, we 
arrive at our conclusions when the probability for a K3 
civilization to expand is of order 10 1 , but they would 
require revision at 10~ 3 . By “expand rapidly” we refer to 
an intergalactic wave of colonization that expands spher¬ 
ically outward from the home galaxy at a substantial 
fraction of the speed of light (we consider here v > .lc), 
generating an expanding cluster of K3 galaxies (use of 
the word “cluster” in this context does not indicate a 
gravitationally bound system merely a localized collec¬ 
tion). This hypothesis is far from self-evident, but we 
will argue in the next section that it should be a natural 
default assumption for K3 civilizations. 

We use the expansion hypothesis in the following way: 
It has recently been shown, in the context of homoge¬ 
neous cosmology, that aggressively expanding civiliza¬ 
tions can rapidly fill the entire universe with life, even if 
such expanders appear rarely m- The timescale for this 
to happen is controlled by the appearance rate and the 
expansion speed of these civilizations. The rate at which 
the universe fills with advanced life in turn influences the 
time of arrival distribution for naturally appearing life 
such as humanity, as the number of untouched “galaxies 
to arrive in” is cut off. Thus, if we specify a scenario by 
the characteristic speed of the expanders, and specify hu¬ 
manity’s relative time of arrival (as mean, 2cr latecomer, 
etc.) then we have fixed the scale of the appearance rate 
for the expanders. The expansion hypothesis then as¬ 
serts that the rate of appearance for all K3 civilizations 
(including non-expanders) should not be orders of mag¬ 
nitude different from this appearance rate, and we are in 
a position to calculate observational quantities like “the 
expected number of K3 civilizations visible within co¬ 
moving radius R” for that scenario. 

Using this approach, we model 81 cosmological scenar¬ 
ios of the type described in m where expansion speeds 
range from .lc to .9c, humanity is regarded as appearing 
at the mean time of arrival or as a la or 2<r latecomer, 
and three different appearance rate functions for intel- 
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ligent life are utilized (depending on the formation rate 
of earthlike planets and assumptions regarding galaxy¬ 
wide extinction events that could delay the appearance 
of advanced life). In nearly all scenarios, we find that 
the co-moving observation distance required to see (on 
average) a single K3 civilization is cosmological, i.e. at 
least as far as the universe’s homogeneity scale of ss .25 
Gly, and much farther in many cases. There also exist 
a significant number of scenarios in which the average 
number of visible civilizations is less than one, no matter 
how far we are able to look - this happens, for example, 
in all scenarios we consider where civilizations expand 
at .9c. The few scenarios that violate this trend, sug¬ 
gesting that many observable K3 civilizations should be 
found nearby, correspond to cases in which humanity has 
arrived improbably late and the expansion speed seems 
improbably slow for a K3-capable civilization. 

When observability is limited to a few Gly (due to the 
practical limitations of a survey), we find that the proba¬ 
bility that K3’s are within visible range is dominated by 
cases in which expanders have entered the visible range 
from farther away. Assumptions regarding galaxywide 
catastrophes (due to gamma ray bursts, etc.) and their 
effects on the appearance rate of life, even unrealistically 
severe ones, have a modest effect on our conclusions. 

This paper is organized in the following way: Section 
II is a brief argument for the expansion hypothesis on 
grounds that are independent from the main analysis of 
this paper. Section III is a review of aggressive expan¬ 
sion scenarios, in which the saturation of the universe by 
advanced life resembles a first-order cosmological phase 
transition involving spatially random “nucleation events” 
followed by spherical expansion. The model of observ¬ 
ability we use is also developed in this section. Section 
IV develops the three life appearance rate models we use, 
which constitute one of the basic inputs of our analysis. 
Section V tabulates the results of our 81 scenarios, orga¬ 
nized by the assumed relative appearance time of human¬ 
ity, while section VI contains a discussion of the results 
and our conclusions. 


II. AN ARGUMENT FOR THE EXPANSION 
HYPOTHESIS 

Although discussions about the possible behavior of 
advanced life tend to be crippled by a severe lack of data, 
in the present context we have the advantage that the ex¬ 
pansion hypothesis refers specifically to K3 civilizations, 
and that carries a number of starting assumptions and 
implications to work with. In particular, the following 
assumptions seem safe: 

1. They have mastered interstellar travel. 

2. They are not adverse to large-scale expansion for 
some fundamental reason. 

3. They place some value on utilizing natural re¬ 
sources at great distances from their origin. 


A recent analysis has made a strong case that inter- 
galactic travel is essentially no more difficult or expen¬ 
sive than interstellar travel - it merely takes longer [2], 
It has also been pointed out that high-speed space travel 
is likely to be the least of the technological hurtles on 
the path to K3 capability, when one considers the re¬ 
quirements implicit in the engineering of entire solar sys¬ 
tems [5]. Our first assumption thus seems to imply that 
practical intergalactic travel should easily be available to 
any K3 civilization. If they have achieved K3 status, then 
they have the means to continuously expand. 

The second assumption also seems to generalize imme¬ 
diately from the case of ss 10 11 stars (a single galaxy) to 
intergalactic travel. Many possible reasons have been 
proposed on the subject of why an advanced civiliza¬ 
tion might choose to stay close to their homeworld and 
focus inward rather than outward [5} H2|, but K3 civ¬ 
ilizations, by definition, have found none of them to 
be compelling. If they have achieved K3 status, they 
cannot be fundamentally inhibited where large-scale ex¬ 
pansion is concerned - they must have expanded ex¬ 
ponentially already on the galactic scale m- Further¬ 
more, a K3 civilization has some motive to utilize re¬ 
sources on a grand scale, following assumption number 3 
- they are not merely neutral on the issue of expansion. 
If such maximally-advanced civilizations have developed 
self-replicating spacecraft so that the cost of such a ven¬ 
ture is minimal, even the mildest preference for expansion 
occurring at any one of the 10 11 solar systems is all that 
will be required. 

Our argument for the expansion hypothesis is essen¬ 
tially that K3 civilizations have, by definition, already 
exhibited all of the technological capability and behavior 
requirements of an aggressive expander, and in the ab¬ 
sence of some powerful, universal, and not-yet-articulated 
reason to stop (or slow dramatically) at the boundaries 
of a home galaxy, it could be assumed that a significant 
fraction of K3’s will continue their expansion at inter¬ 
galactic distances unless constrained by their encounters 
with other expanding civilizations. 


III. AGGRESSIVE EXPANSION SCENARIOS 

An “aggressive expansion scenario” is a proposed cos¬ 
mological phenomenon m, whereby a subset of ad¬ 
vanced life appears at random throughout the universe 
and expands in all directions, saturating galaxies and 
utilizing resources as they go. Mathematically, the de¬ 
scription is almost identical to bubble nucleation and 
growth in a first-order cosmological phase transition, due 
to the common elements of spatially random nucleation 
and spherical expansion. Depending on where the prac¬ 
tical limits of technology lie (in particular, if life is able 
to accelerate the conversion of mass in the universe into 
radiation), heat may also be rapidly released in such a 
scenario, inducing a backreaction on the cosmic scale fac¬ 
tor and pushing the phase transition analogy closer still. 
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Here, we consider a simplified scenario which does not 
take into account heat generation or cosmological back- 
reaction. We also assume that all aggressive expanders 
will be of the same behavior type, i.e. they all expand 
with the same velocity v in the local comoving frame, 
and the expanding spherical front of galaxy colonization 
leads to observable changes a fixed time T after the front 
has passed by. An approximate uniformity of behavior 
of this kind would be expected if the limits of practical 
technology induce an attractor state in the development 
of aggressive expanders throughout the universe, though 
this is only one possibility. 

In such a scenario with uniform expansion behav¬ 
ior, the fraction of the universe that remains unsat¬ 
urated with life, git), can be expressed in Guth-Tye- 
Weinberg [9, 10; (GTW) form as: 

g{t) = e-S~ (i) 


where f(t) is the appearance rate of expanding civiliza¬ 
tions per unit co-moving volume, per unit time, and 
V{t',t) is the volume of space fully saturated with life 
at time f by a single civilization that began expanding 
at t'. The next section will focus on models for f(t). 
When there is a time delay T between the initial arrival 
of expanding spacecraft at some point in space and the 
full saturation of matter there (resulting in observable 
changes), the volume function is given by: 
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where 9(t) is the Heaviside step function and a(t) is the 
cosmic scale factoiQ Given a background cosmological 
solution, then, an aggressive expansion scenario is speci¬ 
fied by giving {v,T, f{t)}. For the scenarios we examine, 
T will play a very minor role in the quantities we calculate 
and could be set to zero as an additional simplification, 
but for the sake of completeness we will take T to cor¬ 
respond to an ample galaxy colonization time [12] of .01 
Gyr, and leave it constant throughout our analysis. 

We now come to a key point of our analysis. Consider 
the set of all human-stage civilizations who will ever have 
appeared in the universe, and who, like humanity, have 
appeared within a non-K3 galaxy. Our prior assumption 
is that we (humanity) are “typical” within this set, and 
in particular that our time of arrival is typical within this 
set. We also assume that the appearance rate for this set 
has the same baseline cosmic time dependence as f{t), 
though the overall proportionality constant could be dif¬ 
ferent by a large (and unknown) constant factor. The 


1 a(t) is taken to be a flat FRW solution with il \o = -683, Q r(l = 
3 X 10 5 , flmO = 1 — Oao — f2 r o, and Ho = .069 Gyr -1 , fixing the 
present age of the universe at to a: 13.75 Gyr. We work in co¬ 
moving coordinates, and use units of Gyr and Gly for dimensions 
of time and distance. 


time of arrival distribution for this set, however, must be 
proportional to the product <?(t)/(f), to account for the 
universe filling up with K3 galaxies, in which no members 
of humanitys set may thereafter appear. In other words, 
the factor of g(t) changes the distribution of arrival times 
significantly by cutting off the appearance rate at late 
times, as the universe fills with advanced life and the pos¬ 
sibility of evolving in an empty galaxy abruptly comes to 
an end. If g(t) = 1 forever (i.e. if there are no aggressive 
expanders), then the distribution and our relative time 
of arrival are fixed by the assumptions going into the 
construction of f(t). Including the possibility of aggres¬ 
sive expansion, however, allows more possibilities and in 
particular we can find (by numerical search) expansion 
scenario parameters that put humanity at to = 13.75Gyr 
at the mean time of arrival or as a 1 or 2 latecomer. The 
mean time of arrival, g, and standard deviation, a, are 
given by: 


fi = N tg{t)f(t)dt 



poo 

/ t 2 g(t)f(t)dt 
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where N = (/ 0 °° g(t)f(t) dt ) 1 for normalization. This 
normalization factor is the reason the unknown propor¬ 
tionality constant for the appearance rate for human- 
stage life does not affect the time of arrival distribution 
one can see that multiplying /(f) by any constant factor 
will not change the mean or standard deviation of arrival 
times. 

The most powerful search to date [8] (by several orders 
of magnitude) for K3 civilizations has involved data from 
full-sky surveys of limited range, so we will be interested 
in calculating EV{obs), the average number (expected 
value) of civilizations that are observable out to some co¬ 
moving distance R. Here, R is supposed to represents a 
limit to one’s equipment and observation techniques (it 
can also be interpreted as the time t\ appearing in fig¬ 
ure 1, and is connected to to through R = dt). 

Regions within the past light cone that can produce an 
observable expanding civilization are illustrated in figure 
1: A corresponds to a region in which any produced K3’s 
are directly observable, while B represents a region in 
which a civilization could appear and, if expanding ag¬ 
gressively, would arrive within A and saturate galaxies 
there, making that civilization observable within R. The 
region C (bounded by a “past saturation cone”) is ex¬ 
cluded because any expanders appearing there would by 
now have fully saturated our own galaxy with advanced 
life, and this is assumed to be ruled out by observation. 
Expanders can thus be seen if they appear in AuB, while 
non-expanding K3’s can be seen if they appear in A U C. 

Considering first only the expanding civilizations, an 
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FIG. 1: The past light cone of an observer at t 0 with 
regions highlighted in which the appearance of an 
expanding civilization is detectable, under the 
assumption that only galaxies out to distance P can be 
directly analyzed for the presence of K3’s. Expanding 
civilizations appearing in region B are observable 
because they expand into region A. Region C is 
presumed to be devoid of expanders because appearing 
in C would imply that our own galaxy is already fully 
saturated with maximally advanced life. 


upper bound for EV ( obs) can be expressed by: 


EV(obs) < 



f(t)V R (t,t 0 )dt 


(5) 


where the volume 


V R (t, t 0 ) of region A U B at time t is: 


Vnit, to) 


0 

4tt ( ft o 1 _ 4 tt 
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This should be regarded as an upper bound be¬ 
cause it includes the possibility of “virtual civilizations” 
- expanding civilizations that appear within already- 
saturated space - and these should not be counted as 
independently observable events. A lower bound on 
EV (obs) can be expressed as: 

EV(obs) > g(t) f(t) V R (t, t 0 ) dt. (7) 

Jo 

This is a lower bound because the fraction of unsaturated 
space in A and B is likely to be higher than g(t), since 
we have additional knowledge that no expanders from C 
have saturated any of the space in A and B. 

For any given expansion scenario (which we will choose 
by fixing the relative arrival time of humanity, and the 
velocity of the expanders), we will want to find a char¬ 
acteristic distance R that represents the power of obser¬ 
vation required to see an expanding civilization. We will 
do this by performing a numerical search of R such that 
EV(obs) = 1, referring to the solution as R±. Because 
EV(obs) = 1, the error contributed by virtual civiliza¬ 
tions will tend to be small, and we will thus use the 
upper bound given above as our estimate for EV(obs) 
in the numerical search of section V. In this approxi¬ 
mation, Ri can also be interpreted as the observation 
distance at which the probability to see zero expanding 


: t > to — T 

(/ t t0 ' T 4r di ') 3 3 : t 0 -T>t> tl (6) 

f(ft°~ T WJ dt ') 3 : h>t 


civilizations is e 1 ~ 37% (the assumption being that 
the appearance of K3’s is a Poisson process, so in general 
P(0) = e ~EV(obs)y 

Although full-sky surveys have been most powerful so 
far, deep field surveys might also be used, which would 
correspond to a small angle in the sky, but an unlimited 
P, i.e. a t\ which extends back to the time when f(t) 
first becomes significantly different from zero. Thus, for 
each of the scenarios we will examine, we also estimate 
EV (obs) (using the upper bound) for the case of an un¬ 
limited R. 

Finally, we can add non-expanding K3 civilizations to 
the analysis. In keeping with the expansion hypothe¬ 
sis, for each scenario considered we will also calculate 
the expected number of visible but non-expanding K3 
galaxies out to Pi under the assumption that the ap¬ 
pearance rate for non-expanders is identical to the ap¬ 
pearance rate for expanders. This corresponds to eval¬ 
uating f*° f(t) V R (t,to) dt where V R (t,to) is the volume 
associated with region AuC from figure 1: 


t 0 ) 


0 


f (L+tAvV) 
^P 3 


3 


t > to — T 
to ~ T > t > *i(8) 
ti > t. 
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IV. APPEARANCE RATE MODELS 


A basic input of an aggressive expansion scenario is the 
appearance rate of expanders per unit coordinate volume, 
per unit cosmic time, f(t). We will consider three such 
models for the time-dependence of f(t), leaving the over¬ 
all proportionality constant as a parameter to be fixed 
by assumptions on the relative time of arrival of human¬ 
ity in the next section. The baseline “non-catastrophic” 
model will set the appearance rate at time t to be propor¬ 
tional to the number of earthlike planets formed between 
4.5 Gyr — 6 Gyr prior to t. This means that we assume 
it takes at least 4.5 Gyr for maximally advanced life to 
appear on a newly formed earthlike planet, and that the 
window for life to evolve is no more than 6 Gyr. This 
assumption is heavily influenced by the successful evo¬ 
lution of intelligence on the Earth, guiding our intuition 
that conditions should be earthlike, and this assumption 
could easily be modified if conditions need not be too 
similar to that of the Earth. The effect of extending the 
“closing of the window” for the evolution of intelligent life 
is to move back the maximum value of f(t), but would 
have little effect on our analysis up to the present time, 
t 0 . Adjusting the time until the “opening of the win¬ 
dow,” however, will shift the initial rise of f{t) in cosmic 
time, though we expect that the opening of the window 
is less likely to vary substantially from our estimate. 

The baseline, non-catastrophic model can be expressed 
as: 

rt- 4.5 

f(t) =a PFR(t') dt' (9) 

Jt -6 


where the planet formation rate, PFR(t ), is modeled by 
PFR(t) = NM(t)SFR(t) with M(t) representing an 
average universe metallicity and SFR(t) the star forma¬ 
tion rate of the universe and N a normalization constant. 
The overall proportionality constant, a, is a free param¬ 
eter to be fixed by time-of-arrival considerations in the 
next section. The buildup of metallicity in the universe 
is, in turn, modeled as an integral over the star forma¬ 
tion rate, M(t) = f* SFR(t') dt\ and we express the star 
formation rate as 


SFR.(t) = 


|lO t_3 : t < 3 

(t-3) 

10-13-75=3 : t > 3. 


( 10 ) 


Here, t is in units of Gyr, representing a simple approx¬ 
imation to the SFR data in m The overall normaliza¬ 
tion for SFR(t) and PFR{t) are chosen such that their 
maximum values are equal to unity. The choice of nor¬ 
malization is essentially arbitrary in this model, corre¬ 
sponding to a rescaling of a. 

This model of PFR(t) is a simplified version of 
Lineweaver’s model m, which additionally takes into 
account a stellar distribution over metallicity. The re¬ 
sulting PFR(t), plotted in figure 2, mirrors the major 
features of that model. 

In addition to the baseline model, we introduce two 
models that include galaxywide extinction events, with 



FIG. 2: Star formation rate SFR(t ), planet formation 
rate PFR(t), and the appearance rate for aggressively 
expanding life f{t) (for a = 1). SFR and PFR are 
normalized to a maximum value of unity. 


a rate that changes as a function of cosmic time. These 
models assume that a life-harboring planet will be sub¬ 
ject to a high-energy event, nearby gamma ray bursts 
(GRB’s) being the prototype example [T|, that destroys 
planetary ozone layers, causing a mass extinction event 
which sets back the evolution of life by some amount of 
time. We will assume that such events are severe if they 
occur in the final stages of evolution towards intelligence, 
so that we reduce the pool of potential planets to those 
which have not seen such an extinction event in the last 
.2 Gyr. This number was chosen to agree with the esti¬ 
mate given by Annis |T], though it may be in the severe 
range of what is plausible. The extinction events are 
modeled as an inhomogeneous Poisson process with in¬ 
tensity A(t) such that A(to) = - this value is chosen 

to be in agreement with a recent analysis jT8j that found 
the probability of a “biospherically important event” oc¬ 
curring on the Earth to be 50% in the last .5 Gyr. 

The two catastrophic models we present differ in their 
assumptions regarding the time-dependence of A (t) - one 
tracks the observed rate of GRB’s in the universe as a 
function of cosmic time, while the other is an extreme 
scenario, with past extinction events far more common 
than suggested by the GRB model. While we expect the 
GRB-tracking model to be realistic, the extreme catas¬ 
trophic model is included to illustrate the extent to which 
our conclusions change under drastic changes to the life 
appearance model, and is not intended to be realistic. 

To be clear, we are not attempting to model all extinc¬ 
tion events in this analysis - only those extinction events 
whose rate we expect to change strongly as a function 
of cosmic time. The rate and effects of local, planetary- 
scale events are assumed to be approximately equivalent 
in cosmic time for sufficiently earthlike conditions, and 
are thus absorbed into the overall proportionality con¬ 
stant a. 

Modeling these catastrophes amounts to multiplying 
the baseline appearance rate model for f(t) with the 
probability that a potential planet has not experienced 
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FIG. 3: The three appearance rate models for /(f) (for 
a = 1), corresponding to the non-catastrophic baseline 
model, the GRB-tracking catastrophic model, and the 
extreme catastrophic model. 


a catastrophe in the last .2 Gyr. For an inhomogeneous 
Poisson process with intensity A(t), this probability is 
e~ /t-. 2 Gyr ) dt The GRB-tracking catastrophic model 
is given by: 


A GRB{t) = £^a(*) 21 (11) 

where the time-dependence comes from the GRB rate 
proportional to (1 + z) 21 found by Wanderman and Pi- 
ran [20]. The extreme catastrophic model is given by: 

A extreme(t) = 6 ~4(*-*o) (12) 

Gyr 

and the time-dependence is chosen arbitrarily. The ef¬ 
fects of these models on /(f) are plotted in figure 3. 

As mentioned in the previous section, if there are no 
expanders, then the time-dependence of /(f) alone (and 
not the proportionality constant a) specifies the time of 
arrival (TOA) distribution of intelligent life. Table I lists 
the mean time of arrival, y, and standard deviation, er, 
for each model. In all three models, the present time, 
to ~ 13.75 Gyr, is a slightly early but completely unsur¬ 
prising time of arrival - this picture will change substan¬ 
tially in the next section when the effect of expanders is 
included. In particular, y and a will be heavily influenced 


by the appearance and expansion rates of the aggressive 
expanders. 



y (Mean) 

a (S.D.) 

Non-Catastrophic 

14.9 Gyr 

5.3 Gyr 

GRB-Tracking Catastrophic 

15.6 Gyr 

5.5 Gyr 

Extreme Catastrophic 

17.1 Gyr 

5.2 Gyr 


TABLE I: Mean times of arrival and standard deviation 
for the three life appearance models, when the effect of 
aggressive expansion is discounted. In all models, the 
appearance time of humanity at to = 13.75 Gyr is typical 
but slightly early. 

V. MODEL RESULTS 

Having described our modeling assumptions and tech¬ 
niques, we are now in a position to numerically examine 
a range of scenarios, looking for their observational con¬ 
sequence. We will divide this into three subsections - 
one for each assumed time of appearance for humanity, 
i.e. to = y, to = y + a, or to = y + 2cr. To reiterate, 
to remains equal to 13.75 Gyr in all cases — it is only 
the relative time of arrival that changes between scenar¬ 
ios (due to the scenario-dependent TOA distribution). 
Then, for each appearance time, we will examine three 
sets of scenarios corresponding to each appearance model 
(non-catastrophic, GRB-tracking catastrophic, and ex¬ 
treme catastrophic). Each set then consists of nine sce¬ 
narios, corresponding to expansion speeds from .lc to 
.9c. 

Each scenario is obtained through a numerical search 
of a required to satisfy the time of arrival condition. Af¬ 
ter each scenario is obtained, we numerically find and 
list the values for R\ such that the expected number of 
observable expanding civilizations is equal to unity, as 
described in section III (also corresponding to the dis¬ 
tance one would have to look to have a probability of 
e^ 1 ~ 37% of seeing zero expanders). The expected num¬ 
ber of visible civilizations, EV(obs ), for unlimited R is 
also reported, as will be the expected value of observable 
non-expanders (isolated K3 galaxies), under the assump¬ 
tion that their appearance rate is identical to that of the 
expanders (in keeping with the expansion hypothesis). 

In a significant number of cases, particularly associated 
with the mean time of arrival condition, EV ( obs ) is less 
than unity for unlimited R (i.e. R\ is undefined), due 
to the expanders having being made extremely rare to 
satisfy the time of arrival condition. These cases will be 
marked with N/A for the relevant quantities. 
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A. Humanity (to = 13.75 Gyr) at the mean time of arrival 

1. Non-catastrophic appearance model 


Expansion Speed 

v = .lc 

v = .2c 

v = .3c 

v — Ac 

v = .5c 

v = .6c 

v = .7c 

v = .8c 

v = .9c 

a (appearances/Gly 3 Gyr) 

.019 

.0024 

.00071 

.00030 

.00015 

.000089 

.000056 

.000038 

.000026 

Ri (Gly) 

1.0 

2.4 

5.0 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

EV ( obs ) for unlimited R 

35 

4.3 

1.2 

.51 

.24 

.13 

.066 

.033 

.013 

EV (obs) of non-expanders within Ri 

.52 

.69 

.95 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 


TABLE II: Mean TOA, non-catastrophic appearance model 

2. GRB-tracking catastrophic appearance model 


Expansion Speed 

v = .lc 

v = .2c 

v = .3c 

v = Ac 

v = .5c 

v = .6c 

v = .7c 

v = .8c 

v = .9c 

a (appearances/Gly 3 Gyr) 

.055 

.0069 

.0020 

.00086 

.00044 

.00025 

.00016 

.00011 

.000075 

Ri (Gly) 

.77 

1.8 

3.4 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

EV(obs) for unlimited R 

57 

7.1 

2.1 

.84 

.40 

.21 

.11 

.054 

.021 

EV (obs) of non-expanders within Ri 

.44 

.58 

.79 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 


TABLE III: Mean TOA, GRB-tracking catastrophic appearance model 

3. Extreme catastrophic appearance model 


Expansion Speed 

v = .lc 

v = .2c 

v = .3c 

v = .4c 

v = .5c 

v = .6c 

v = .7c 

v = .8c 

v = .9c 

a (appearances/Gly 3 Gyr) 

.54 

.068 

.020 

.0085 

.0043 

.0025 

.0016 

.0011 

.00075 

Ri (Gly) 

.40 

.92 

1.7 

2.9 

N/A 

N/A 

N/A 

N/A 

N/A 

EV(obs) for unlimited R 

110 

13 

3.9 

1.6 

.75 

.39 

.21 

.10 

.040 

EV (obs) of non-expanders within R,\ 

.34 

.45 

.61 

.87 

N/A 

N/A 

N/A 

N/A 

N/A 


TABLE IV: Mean TOA, extreme catastrophic appearance model 






B. Humanity (to = 13.75 Gyr) as a la latecomer 


1. Non-catastrophic appearance model 


Expansion Speed 

v = .1 c 

v = .2c 

v = .3c 

v = Ac 

v = .5c 

v = .6c 

v = .7c 

v = .8c 

v = .9c 

a (appearances/Gly 3 Gyr) 

.25 

.031 

.0093 

.0039 

.0020 

.0012 

.00073 

.00049 

.00034 

Ri (Gly) 

.24 

.55 

.94 

1.5 

2.3 

3.6 

N/A 

N/A 

N/A 

EV(obs) for unlimited R 

450 

56 

16 

6.6 

3.2 

1.6 

.86 

.43 

.17 

EV(obs) of non-expanders within Ri 

.10 

.14 

.20 

.30 

.49 

.85 

N/A 

N/A 

N/A 


TABLE V: a latecomer, non-catastrophic appearance model 


2. GRB-tracking catastrophic appearance model 


Expansion Speed 

v = .lc 

v = .2c 

cS 

II 

CO 

o 

c2 

II 

v = .5c 

II 

o 

c2 

II 

O 

°o 

II 

v = .9c 

a (appearances/Gly 3 Gyr) 

.51 

.063 

.019 

.0079 

.0040 

.0023 

.0015 

.00099 

.00069 

Ri (Gly) 

.21 

.47 

.82 

1.3 

2.0 

3.1 

7.3 

N/A 

N/A 

EV(obs) for unlimited R 

520 

65 

19 

7.7 

3.7 

1.9 

1.0 

.50 

.19 

EV (obs) of non-expanders within Ri 

.090 

.12 

.18 

.27 

.43 

.74 

1.5 

N/A 

N/A 


TABLE VI: a latecomer, GRB-tracking catastrophic appearance model 


3. Extreme catastrophic appearance model 


Expansion Speed 

v = .lc 

v = .2c 

v = .3c 

v = .4c 

v = .5c 

v = .6c 

v = .7c 

v = .8c 

v = .9c 

a (appearances/Gly 3 Gyr) 

3.7 

.46 

.14 

.057 

.029 

.017 

.011 

.0072 

.0050 

Ri (Gly) 

.12 

.27 

.47 

.73 

1.1 

1.7 

3.0 

N/A 

N/A 

EV(obs) for unlimited R 

720 

90 

26 

11 

5.1 

2.6 

1.4 

.69 

.27 

EV (obs) of non-expanders within Ri 

.070 

.096 

.14 

.21 

.33 

.57 

1.1 

N/A 

N/A 


TABLE VII: a latecomer, extreme catastrophic appearance model 
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C. Humanity (to = 13.75 Gyr) as a 2 a latecomer 

1. Non-catastrophic appearance model 


Expansion Speed 

v — .lc 

v = .2c 

v = .3c 

v = Ac 

v = .5c 

v = .6c 

ci 

II 

O 

°o 

II 

v = .9c 

a (appearances/Gly 3 Gyr) 

.92 

.12 

.034 

.014 

.0074 

.0043 

.0027 

.0018 

.0013 

Ri (Gly) 

.089 

.20 

.35 

.54 

.81 

1.2 

1.9 

3.5 

N/A 

EV(obs) for unlimited R 

1700 

210 

60 

24 

12 

6.0 

3.2 

1.6 

.61 

EV(obs ) of non-expanders within Ri 

.019 

.027 

.040 

.062 

.10 

.19 

.42 

1.2 

N/A 


TABLE VIII: 2 a latecomer, non-catastrophic appearance model 

2. GRB-tracking catastrophic appearance model 


Expansion Speed 

v = .lc 

v = .2c 

v = .3c 

CJ 

II 

v = .5c 

v = .6c 

v = .7c 

<5 

II 

oo 

v = .9c 

a (appearances/Gly 3 Gyr) 

1.8 

.22 

.066 

.028 

.014 

.0082 

.0052 

.0035 

.0024 

Ri (Gly) 

.080 

.18 

.30 

.47 

.71 

1.1 

1.7 

3.0 

N/A 

EV(obs) for unlimited R 

1800 

230 

67 

27 

13 

6.7 

3.5 

1.8 

.68 

EV (obs) of non-expanders within Ri 

.017 

.024 

.035 

.055 

.091 

.17 

.37 

1.1 

N/A 


TABLE IX: 2 a latecomer, GRB-tracking catastrophic appearance model 


3. Extreme catastrophic appearance model 


Expansion Speed 

v — .lc 

v = .2c 

v = .3c 

II 

v = .5c 

<5 

II 

Oi 

O 

v = .7c 

v = .8c 

v = .9c 

a (appearances/Gly 3 Gyr) 

13 

1.6 

.48 

.20 

.10 

.060 

.038 

.025 

.018 

Ri (Gly) 

.044 

.098 

.17 

.26 

.39 

.59 

.92 

1.6 

N/A 

EV(obs) for unlimited R 

2600 

320 

92 

37 

18 

9.3 

4.9 

2.4 

.95 

EV (obs) of non-expanders within Ri 

.012 

.016 

.024 

.037 

.063 

.12 

.25 

.73 

N/A 


TABLE X: 2cr latecomer, extreme catastrophic appearance model 


VI. DISCUSSION AND CONCLUSIONS 

From the model results of the previous section, some 
conclusions are immediate: 

• There exist large regions of the parameter space 
that result in EV(obs) < 1 for unlimited R. The 
probability must be regarded as substantial that we 
can see no expanders or K3 civilizations, no matter 
how good our observation techniques, even if the 
saturation of the universe by advanced life is well 
underway. 

• If humanity has appeared near (or before) the 
mean time of arrival for civilizations like ours, the 
prospect of seeing any expanders or isolated K3 civ¬ 
ilizations seems poor. For observation to be likely 
at the mean arrival time, the limits to technology 


would have to make intergalactic expansion prac¬ 
tical, but not above « .3c. Other authors have 
concluded that such a barrier could be surpassed 
even by relatively simple fission rockets [2], mak¬ 
ing such low-u scenarios seem less plausible. If we 
are living in such a scenario and a positive obser¬ 
vation is made, it would most likely be at distances 
of multiple Gly. 

• If humanity is a ler latecomer, observational 
prospects are better. Expansion scenarios up to 
.7c are likely to be observable, at least in principle. 
For such scenarios, R\ remains cosmological and at 
multiple Gly when expansion above .4c is practical 
for maximally advanced life. 

• In the seemingly unlikely case that humanity is a 
2cr latecomer, prospects for observation are good in 
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expansion scenarios up to nearly .9c. In fact, the 
low-c scenarios of this case are probably ruled out 
already by existing observations {5j, provided they 
can be regarded as sufficiently thorough searches 
for K3 civilizations. Even if we are a 2a latecomer, 
we still expect to make the first observations at 
cosmological distances greater than the homogene¬ 
ity scale (if expansion speeds are greater than .2c), 
though only expansion speeds above .6c result in 
an Ri above 1 Gly. 

• In no cases examined do we find expansion scenar¬ 
ios at .9c such that R\ is defined. However, if we are 
a 2er latecomer in the GRB-tracking catastrophic 
scenario, the probability that zero expanders are 
visible is e ~ EV ( obs ) « 51%. We should not realis¬ 
tically expect to see K3’s if intergalactic travel at 
v = .9c is practical and the expansion hypothesis 
is correct, but the probability is not overwhelm¬ 
ingly negative. The difficulty in observing high-u 
scenarios comes from two factors - the very small 
appearance rate required to satisfy the time of ar¬ 
rival conditions when the expansion speed is high, 
and the fact that in high-u scenarios, a large frac¬ 
tion of our past light cone (region C of fig. 1) is 
already known to be devoid of such expanders (else 
they would already be here). 

• The expected number of non-expanders appearing 
within Ri is almost always less than the number of 
expanders that are visible out to R± (unity), despite 
the fact that non-expanders had the additional op¬ 
portunity to appear in region C of fig. 1. Due to 
the assumed equality of appearance rates (between 
expanders and non-expanders), this can be inter¬ 
preted to mean that we are more likely to observe 
an expander that came to within Ri from farther 
away, rather than seeing an expander that appeared 
within Ri to begin with. It also means that we ex¬ 
pect to have more K3 clusters within visible range 
than isolated K3 galaxies. 

• Even unrealistically extreme models of galaxy-scale 
extinction events have a modest effect on our con¬ 
clusions. The main effect of such models is to make 
advanced life nearly impossible in the early uni¬ 
verse, before advanced life would have time to arise 
anyway. The effect is diminished in more recent 
times, which are more relevant for the appearance 
of advanced life. 

It is interesting to visualize just how rarely aggressively 
expanding civilizations arise, according to this analysis. 
A typical value for the appearance rate parameter a in a 
GRB-tracking scenario is of order ICR 3 appearances per 
Gly 3 per Gyr. In other words, it would take a sphere of 
radius ~ 5 Gly to produce a single aggressive expander 
in a billion years. This is a volume encompassing many 


thousands of superclusters and perhaps a hundred million 
large galaxies. Similar numbers have been implied by [2, 
[7], in their calculation of the number of galaxies that 
could have reached and colonized the Milky Way. The 
great filter im implied by this type of universe must be 
very great indeed. 

We should reflect on our use of the time of arrival dis¬ 
tribution as a means of estimating the appearance rate 
parameter a. This is a form of anthropic reasoning, im¬ 
plicitly utilizing the Self-Sampling Assumption, which 
exhorts us to reason as if we are a random sample from 
the set of all comparable observers who will ever have 
existed 0. We have interpreted this to mean that our 
time of arrival should be typical in the set of human- 
stage civilizations who will have appeared in a non-K3 
galaxy. Though we have no prior theory to determine a, 
we should feel confident that a scenario in which human¬ 
ity is a multiple-cr latecomer to the universe can be re¬ 
garded as very improbable without good evidence to the 
contrary. This allows us to focus our attention on sce¬ 
narios where humanity is no more than a 2a latecomer. 

Limitations of our analysis should also be noted. Since 
we have assumed a homogeneous universe, the visible ge¬ 
ometry of small domains (smaller than the homogeneity 
scale of the universe) could deviate significantly from the 
expanding spheres assumption. Similarly, the universe- 
averaged appearance rate is a rough approximation, at 
the current level of development. Taking into account the 
details of galaxy evolution could presumably make signif¬ 
icant changes to /(f), though we have seen that our con¬ 
clusions seem to remain fairly robust with respect to sub¬ 
stantial changes to /(f). The assumption of behavior uni¬ 
formity (i.e. a single, constant {v,T} for the expanders) 
is also debatable - a mixture of rare fast expanders and 
common slow expanders, for example, might be expected 
to change our conclusions substantially, based on our 
previous calculations m as would extreme galaxy col¬ 
onization models that take T to be on the scale of mul¬ 
tiple Gyr. We do not expect the remaining uncertainty 
in the underlying background cosmological parameters 
to significantly affect our conclusions, provided that the 
A CDM model remains standard. 

From the most practical point of view, what does our 
analysis say? Our modeling is consistent with the possi¬ 
bility that K3 civilizations and aggressive expanders are 
present but not observable, or even not present at all. 
But we do have a conditional result - if K3 galaxies 
are observable, then we expect to see them at cosmolog¬ 
ical distances, as part of a cluster of K3 galaxies whose 
boundary expands at a middling fraction of the speed of 
light. Conditional results of this kind have a curious sort 
of practicality: Knowing the location of a streetlight does 
not tell us where we have dropped our keys, but it does 
let us know where we have the best chance of finding 
them. In this way, such an analysis can be of practical 
use for future searches of K3 galaxies. 


11 


[1] J Annis. An astrophysical explanation for the great 
silence. Journal of the British Interplanetary Society, 
52:19-22, 1999. 

[2] Stuart Armstrong and Anders Sandberg. Eternity in 
six hours: Intergalactic spreading of intelligent life and 
sharpening the fermi paradox. Acta Astronautica, 89:1 
13, 2013. 

[3] Nick Bostrom et al. Anthropic bias: Observation selection 
effects in science and philosophy. Routledge, 2002. 

[4] Robert J Bradbury, Milan M Cirkovic, and Geroge 
Dvorsky. Dysonian approach to seti: a fruitful middle 
ground? Journal of the British Interplanetary Society, 
64:156-165, 2011. 

[5] Richard A Carrigan Jr. Is interstellar archeology possi¬ 
ble? Acta Astronautica, 78:121-126, 2012. 

[6] MM Cirkovic. Against the empire. Journal of the British 
Interplanetary Society, 61:246-254, 2008. 

[7] Martyn Fogg. Feasibility of intergalactic colonisation and 
its relevance to setl. Journal of the British Interplanetary 
Society , page 491, 1988. 

[8] Roger Griffith, Jason Wright, Jessica Maldonado, 
Matthew S Povich, Steinn Sigurdjsson, and Brendan 
Mullan. The g infrared search for extraterrestrial civiliza¬ 
tions with large energy supplies, iii. the reddest extended 
sources in wise. The Astrophysical Journal Supplement 
Series, 217(2):25, 2015. 

[9] Alan H Guth and S-HH Tye. Phase transitions and mag¬ 
netic monopole production in the very early universe. 
Physical Review Letters, 44(10):631, 1980. 

[10] Alan H Guth and Erick J Weinberg. Cosmological conse¬ 
quences of a first-order phase transition in the s u 5 grand 
unified model. Physical Review D, 23(4):876, 1981. 

[11] Robin Hanson. The great hlter-are we almost past it? 
preprint available at http://hanson. gmu. edu/greatfilter. 
html, 1998. 

[12] Michael H Hart. Explanation for the absence of extrater¬ 
restrials on earth. Quarterly Journal of the Royal Astro¬ 
nomical Society, 16:128, 1975. 

[13] Nikolai S Kardashev. Transmission of information by 
extraterrestrial civilizations. Soviet Astronomy, 8:217, 


1964. 

[14] Brian C Lacki. Seti at planck energy: When parti¬ 
cle physicists become cosmic engineers. arXiv preprint 
arXiv:1503.01509, 2015. 

[15] Charles H Lineweaver. An estimate of the age distri¬ 
bution of terrestrial planets in the universe: quantifying 
metallicity as a selection effect. Icarus, 151(2):307-313, 
2001 . 

[16] William I Newman and Carl Sagan. Galactic civiliza¬ 
tions: Population dynamics and interstellar diffusion. 
Icarus, 46(3):293-327, 1981. 

[17] S Jay Olson. Homogeneous cosmology with aggressively 
expanding civilizations. Classical and Quantum Gravity, 
32(21):215025, 2015. 

[18] Tsvi Piran and Raul Jimenez. Possible role of gamma ray 
bursts on life extinction in the universe. Physical review 
letters, 113(23):231102, 2014. 

[19] Carl Sagan and William I Newman. The solipsist ap¬ 
proach to extraterrestrial intelligence. Quarterly Journal 
of the Royal Astronomical Society, 24:113, 1983. 

[20] David Wanderman and Tsvi Piran. The luminosity func¬ 
tion and the rate of swift’s gamma-ray bursts. Monthly 
Notices of the Royal Astronomical Society, 406(3):1944- 
1958, 2010. 

[21] JT Wright, RL Griffith, S Sigurdsson, MS Povich, and 
B Mullan. The g infrared search for extraterrestrial 
civilizations with large energy supplies, ii. framework, 
strategy, and first result. The Astrophysical Journal, 
792(1):27, 2014. 

[22] JT Wright, B Mullan, S Sigurdsson, and MS Povich. The 
g infrared search for extraterrestrial civilizations with 
large energy supplies, i. background and justification. 
The Astrophysical Journal, 792(1):26, 2014. 

[23] Erik Zackrisson, Per Calissendorff, Saghar Asadi, and 
Anders Nyholm. Extragalactic seti: The tully-fisher re¬ 
lation as a probe of dysonian astroengineering in disk 
galaxies. The Astrophysical Journal, 810(1):23, 2015. 



